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ABSTRACT 

In a flat universe dominated by dark energy, the Integrated Sachs- Wolfe (ISW) effect 
can be detected as a large-angle cross-correlation between the CMB and a tracer of 
large scale structure. We investigate whether the inconclusive ISW signal derived from 
2MASS galaxy maps can be improved upon by including photometric redshifts for the 
2MASS galaxies. These redshifts are derived by matching the 2MASS data with optical 
catalogues generated from SuperCOSMOS scans of major photographic sky surveys. 
We find no significant ISW signal in this analysis; an ISW effect of the form expected 
in a ACDM universe is only weakly preferred over no correlation, with a likelihood 
ratio of 1.5:1. We consider ISW detection prospects for future large scale structure 
surveys with fainter magnitude limits and greater survey depth; even with the best 
possible data, the ISW cross-correlation signal would be expected to evade detection 
in > 10% of cases. 

Key viTords: 



1 INTRODUCTION 

It is now widely accepted that the universe today is dom- 
inated by Dark Energy. Evidence for this comes from a 
number of sources: supernovae observations indicate a lu- 
minosity distanc e-redshift relation consistent with acceler - 
ated expansion (|Riess et al.lll998l : [Peilmutter et allll999^ : 
the combination of CMB and galaxy clustering measure- 
ment s places tight constra i nts on the Dark Energy den- 
sity (|Efstathiou et al.lll990l . I2OO2I : ISpergel et alll200'<1 ) and 
Baryon Acoustic Oscillation (BAO) measurements suggest 
an angular diameter distance-redshift rela t ion consistent 
with Dark Energy l|Eisenstein et al.] l2005l : IPercival et al.l 
l200'iD . 

Further independent evidence for Dark Energy can be 
provided by detection of the Integrated Sachs- Wolfe (ISW) 
effect, a gravitational secondary anisotropy in the CMB. The 
ISW effect is a net change in the energy of a CMB photon 
as it passes through evolving gravitational potential wells. 
The ISW temperature fluctuations introduced via the ISW 
effect can be calculated from 



of the photon at time t and $ i s the gravitational potential 
(|Martinez-Gonzalez et al.|[l990l ). 

The ISW effect is a large scale phenomenon: tempera- 
ture fluctuations tend to cancel on smalle r scales as pho- 
tons pass through many potential wells (|Hu fc DodelsonI 
|2002| ). The ISW effect is sub-dominant to the primary CMB 
anisotropics and exists on scales where cosmic variance is 
large, making its detection from CMB data alone infeasi- 
ble. Instead, the cross-correlation bet ween CMB maps and 
trace rs of large scale structure is used (|Crittenden fc Turokl 
ll996|V_ Many different tracers ha ve been correlated e.g. X- 
ray (BouKhn & Crittcndcn"2005'); Sloan Digital Sky Survey 
(Fosalba ct al. 2003; Scranton et al. 2003); Lu minous Red 
Gala pdes (iCranett et al.l l2009l: ISawangwit et all I2OI0I ): ra- 



Tc 



At, 



(1) 



where to and tLs denote the times today and at last scat- 
tering respectively; x is the position along the line of sight 
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dio (iNolta et al. 2004 : Pietrobon et all 20061 ): and quasars 
(jCiannantonio et al.ll20od ). Not all of these lead to signifi- 
cant signals, but the stron gest claimed detect i ons cu rrently 
lie at the 3(j — 4a" level. iGiannantonio et al.l (|2008) cross- 
correlated data from a number of different large scale struc- 
ture surveys with WMAP CMB data and reported an overall 
IS W effect that wa s signi ficant at the 4.5(7 level. 

lAfshordi et al.1 (|2004l ) cross-correlated the Two Micron 
AU S ky Survey (2MASS. |Jarrett|[2004i ) and year-1 WMAP 
data ([Bennett et al.ll2003l ) and claimed to detect the ISW 
effect at 2.5(t. Howe ver, no significant ef fect was found in 
the updated work of iRassat et al.1 (|2007l ) (hereafter R07), 
who correlated 2M AS S with the year-3 WMAP CMB data 
(jHinshaw et al.ll2007h . In the present paper, we investigate 
whether the addition of photometric redshifts to the 2MASS 
survey data can boost the significance of any ISW signal. We 
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also consider detection prospects for similar surveys with 
fainter magnitude limits and surveys with greater depth. 

For ease of comparison with R07 we adopt the 
same cosmological model, namely a flat universe with 
f^m = 0.3, h = 0.7 where the Hubble parameter is 
Ho =iOO h kms~^Mpc~^, spectral index n = 1, Q,i, = 0.05 
and (Tg = 0.75. Section [2] describes the data that we use in 
this paper and details of the masks for the galaxy and CMB 
maps. SectionOoutlines the theoretical cross-correlation sig- 
nal expected in a ACDM universe and Section[4]presents the 
results of our harmonic space cross-correlation analysis. Sec- 
tion [S] assesses the statistical power of the method used here 
as a detection tool, and investigates the detection prospects 
of more ambitious large scale structure surveys. Our conclu- 
sions are presented in Section [S] 



2 DATA 

The large scale structure dataset used in the cross- 
correlation analysis is by far the greatest source of uncer- 
tainty. On the s cales of the ISW effec t, the CMB data are sig- 
nal dominated feperge l et al.l [20071 ) whilst the galaxy data 
are afflicted by shot noise, uncertainties in the bias relation 
and photometric redshift errors. 
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Figure 1. A comparison of photometric and spectroscopic red- 
shifts for 2MASS galaxies. The photometric redshifts are com- 
puted from the BRJHK dataset resulting from combining 2MASS 
with SupcrCOSMOS optical photometry. 



2.1 2MASS galaxy data 

2MASS is an all-sky survey in the J, H and Ks, bands. The 
final extended source catalogue (XSC) contai ns over 1.6 mi l- 
lion objects, over 98% of which are galaxies ( Jarrett| |2004h . 
This dataset is ideal for constructing maps of galaxy mass 
density since near infrared selection means that 2MASS is 
sensitive to old stars and hence the most massive structures. 

Photometric redshifts for the 2MASS XSC have been 
generated by matching the 2MASS data with optical cat- 
alogues generated from SuperCOSMOS scans of the major 
photographic sky surveys (UKST in the south; POSS2 in the 
north). Details of the S uperCOSMOS cata logue construc- 
tion process are given in lHamblv et al. These cata- 
logues have been given an accurate photometric calibration 
using a mixture of SDSS photometry, plate overlaps, and 
requiring uniformity in average colour between optical band 
and 2MASS J ; the details of the process are described in 
iPeacock et al.l (|2010l ). The photographic data define an all- 
sky optical galaxy catalogue in photographic B,R, I bands, 
reaching completeness limits of approximately B = 21.5, 
R = 20.5, / = 19. Since the /-band plates are significantly 
less deep, they were not used in the current analysis, leaving 
a 5-band BRJHK dataset from which photometric redshifts 
were to be derived. Owing largely to SDSS, 2dFGRS, and 
6dFGS, a very substantial set of calibrating spectroscopy 
exists - for about 30% of the catalogue. It was therefore pos- 
sible to take a highly direct approach to the generation of 
photometric redshifts, averaging over neighbours of known 
redshift at a given location in the 5-band magnitude space. 
By using magnitudes rather than colours, this automatically 
builds in information from the luminosity function, so that 
bright galaxies are never allocated an extreme redshift that 
would require them to have an unrealistic luminosity. As a 
result, the scatter in photometric redshift declines towards 
z = 0, as shown in Fig. 1. The overall rms in Zphot — z: is 



0.033. By construction, this process yields photometric red- 
shifts that are unbiased in that the mean true redshift at 
given Zphot is equal to Zphot. An inevitable consequence is 
that there is then a bias in Zphot at given true z: Zphot is 
overestimated near z = Q and underestimated at high z, by 
an offset of 0.02 to 0.03. In any case, what we will need below 
is the ability to predict the distribution of true z from a band 
of Zphot; provided this is known, any bias is unimportant. 

The optical filter properties in UKST and POSS2 are 
very slightly different, so the process was performed sep- 
arately in the two hemispheres. However, there is deeper 
calibrating spectroscopy in the north, which introduces the 
possibility of a small bias. We therefore made a small rescal- 
ing as a function of redshift so that the global redshift dis- 
tributions in each hemisphere were forced to be identical 
(in shape, but not in number) for our standard extinction- 
corrected if-selected galaxy sample. 

We construct 3 galaxy subsamples in approximately in- 
dependent photometric redshift slices of thickness Az = 0.1 
out to 2 = 0.3. The true redshift distribution of galaxies 
in each of these slices is calculated using the calibrating 
spectroscopic information; these are shown in Fig. (2] We 
parametrize the redshift distributions using 

^Kz'^expj-lz/z,)"}, (2) 

and perform a least squares fit to obtain the distribution 
parameters. The best-fitting values of q, /3 and z* for each 
photometric redshift slice are given in Table [1] It is clear 
that the 0.2 < z < 0.3 slice makes very little contribution 
to the 2MASS sky distribution; yet by virtue of its larger 
volume, this slice will double the total ISW signal compared 
to a limit at z < 0.2. We may thus hope that the use of 
redshift information in the current study will improve the 
chances of detection. 
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which models for synchrotron, free-free and dust emission 
have been subtracted. It is however prudent to mask the 
WMAP maps at low Galactic latitudes and the WMAP Kp2 
mask is used to this end, leaving ~ 85% of the sky unmasked. 

The fact that both the galaxy maps and the CMB maps 
are masked will mean that spherical harmonic coefficients 
measured from the maps are not exact and an estimator 
needs to be used. Here, following R07, we compensate for 
the lack of sky coverage using the factor /sky in equation 

m- 

We now summarise the necessary theory of the ISW 
effect (see e.g. Cooray et al. 2002, Afshordi et al. 2004 for 
details). 



Figure 2. The redshift distribution of galaxies selected by pho- 
tometric redshift to lie in slices z = 0.0 — 0.1, z = 0.1 — 0.2 and 
2 = 0.2 — 0.3. The histograms show the distribution of galaxies 
in the photometric redshift shells with spectroscopic information; 
here the redshift distribution functions are normalized to match 
these data. 



Slice 


-^gals 




a 


/3 


z. 


0.0 < z < 0.1 


451329 


195789 


0.88 


2.24 


0.076 


0.1 < z < 0.2 


258673 


101912 


3.46 


2.27 


0.104 


0.2 < z < 0.3 


19162 


8543 


17.4 


2.12 


0.083 



Table 1. The redshift distribution of galaxies in each photomet- 
ric redshift slice is described by equation ([2)l. The table gives the 
total number of galaxies in each slice, the number having spec- 
troscopic redshifts, and best-fitting values of a, fi and z» for each 
photometric redshift slice. 



It is necessary to mask the region of sky around the 
Galactic plane before making the galaxy overdensity maps, 
since the survey is incomplete here and con fusion is more 
likely. Following R07, we use the dust maps of lSchlegel et al.l 
l|l998h to mask regions with K-hanA reddening Ak > 0.05, 
which leaves approximately 67% of the sky unmasked. For 
the cross-correlation analysis we impose magnitude cuts 
12.0 < K < 13.8 on the extinction corrected 2MASS data 
to ensure uniform coverage across the whole sky. The num- 
ber of galaxies remaining in each of the photometric redshift 
slices is also given in Table [T] 

We use HEALPix software (Gorski et al. 2005; see also 
http://healpix.jpl.nasa.gov) to generate maps of the 
galaxy overdensity and to compute the spherical harmonic 
coefficients aim used in the cross-correlation analysis. Fig. 
[3] shows the galaxy overdensity in each of our photometric 
redshift slices. 



2.2 WMAP CMB data 

We take our C MB temperature da ta from the third year 
WMAP resuhs (|Hinshaw et al.ll2007f ). Cross-correlations are 
computed for the 3 bands least contaminated by Galactic 
foregrounds (Q, V and W) and also for the Internal Linear 
Combination (ILC) map. To minimise the effects of Galac- 
tic foregrounds we use the foreground reduced sky maps in 



3 THEORY AND GALAXY BIAS 

Expanding the projected galaxy density field in a given red- 
shift shell and the ISW temperature fluctuations in spherical 
harmonics, the cross correlation between the galaxy overden- 
sity field and the CMB temperature field is; 



(3) 



where af^ and aJm are the spherical harmonic coefficients 
for the galaxy overdensity and temperature fields respec- 
tively and 5^ is the Kronecker delta. 

In Fourier space and within the linear regime, Poisson's 
equation relates the time derivative of the potential field to 
the matter density field; 



dt 



9{a) 



4(^ = 0), 



(4) 



where g{a) is the linear growth factor, normalised to g — 1 
today. From the definition of the ISW temperature fiuctua- 
tions we then have; 



T 



n 

(27r)3c3 



m 

27r2c3 
d^k 



dr a 



1[£] 

dt lal 



— 5kjt{kr)Ylm{^), 



(5) 



where 5 denotes the matter overdensity field today, k and n 
are unit vectors and the exponential has been re-expressed 
in terms of spherical Bessel functions ji and spherical har- 
monics Ylm- 

Within each redshift slice we calculate the projected 
galaxy overdensity (5proj, which is related to the three- 
dimensional galaxy overdensity via 

n(r-)r^ 



jn{r 



dr' 



-5{r) dr= e(r)<5(r) dr. 



(6) 



where n(r) is the background galaxy number density and 
r{z) is comoving distance; thus 0{r) is the effective radial 
weight applied to the density field. 

Assuming a linear bias relation to relate galaxy and 
matter overdensities 5g — bSm, the galaxy expansion coeffi- 
cients are; 
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Figure 3. Maps of the 2MASS galaxy overdensity in each of the 3 redshift shoes considered: 0.0 < z < 0.1 (top), 0.1 < z < 0.2 (middle) 
and 0.2 < z < 0.3 (bottom). The density field is smoothed with a Gaussian of FWHM = 100' for the z = 0.0 - 0.1 and z = 0.1 - 0.2 
slices and of FWHM = 200' for the z = 0.2 - 0.3 sUce. 
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3.1 Galaxy bias 




Figure 4. The expected cross-correlation signal in fiK for each 
of the 3 redshift slices with the Limber approximation (dashed) 
and exact (solid) lines. The peak of the cross-correlation shifts 
to higher multipoles for higher redshift slices (red: < 2 < 0.1; 
green: 0.1 < z < 0.2; blue: 0.2 < z < 0.3). Here a constant bias 
b = 1.4 has been used. 



J drgOir) J 



^ I drgOir) / d^k ^^jf (fer)y^™(k), 



(7) 



where once again the matter overdensity S is evaluated today 
and k is a unit vector. 

Now the cross correlation can be written as 



where / = dlng/dlna. 

Limber's approximation can be used to eliminate the 
spherical Bessel functions by using the small angle limit 



dfc P{k) J dr' ge{r')je{kr') 
X ldrgH{f-l)jt{kr), (8) 



Applying this to equation ((8]) gives 



(9) 



1+1/2 



drQ{r)Hg\f 



xP 



(10) 



Fig. |4] shows the expected cross-correlation for each of 
our galaxy redshift shells. The peak of the cross correla- 
tion can be seen to shift to higher multipoles at higher red- 
shifts as the power spectrum peak shifts to higher multi- 
poles (a given physical scale will subtend a smaller angle 
when viewed from further away, i.e. will correspond to a 
higher multipole). The Limber approximation is valid here 
for £ > 10. 

Following an identical method, the galaxy auto- 
correlation is given by 

|2 



21/ 

TV 



dk k^P{k) 



dr gO{r)ji{kr) 



(11) 



or on applying the Limber approximation l|Kaiserlll992l ). 



dr — P ' 



(12) 



The galaxy bias is expected to change with z: the furthest 
redshift shell is populated only by the brightest, most mas- 
sive galaxies, which are known to be mo re strongly clustered 
(|Park et al.l Il994 iLovedav et"al] 119951 ). In order to deter- 
mine the bias in each redshift slice, we use the galaxy angular 
power spectra. The form of the angular power spectrum in 
each slice can be predicted given a three-dimensional power 
spectrum using equations (|lip and (|12p . The effective bias 
of the predicted spectrum can then be adjusted to match 
observations. In linear theory, bias cannot be determined in- 
dependently of (Tg which also acts to renormalize the power 
spectrum, Ct oc {bag)^. In what follows, we therefore fix the 
value of (Tg at 0.75 an d fit only for 

We use CAMB ijLewis et all l200d ) to generate non- 
linear matter power spectra for our cosmological model (see 
Section [l| and output the three-dimensional spectra at red- 
shifts at the mid-point of each redshift shell in question. 
Although only strictly valid in linear theory, we evolve the 
non-linear power spectrum for each slice according to the 
growth function; this is thought to be a good approxima- 
tion to the non-linear power spectrum on the linear and 
quasi-linear scales where we will fit for the bias. We have 
two sources of data that can be used to measure the an- 
gular power in each redshift slice: the photometric redshift 
data and the subset of this in the SDSS region where com- 
plete spectroscopic redshifts exist. Since the spectroscopic 
data exist only in a subregion, the angular power measured 
in this case will not be accurate on large scales. However, 
we should expect reasonable agreement between the photo- 
metric angular power and the spectroscopic angular power 
on smaller scales, up to differences due to sample variance 
and mask. Fig. [5]shows the shot noise corrected photometric 
and spectroscopic angular power spectra binned for £ ^ 60 in 
bins of width A£ = 10 and the best-fitting predicted spectra 
for each of the redshift slices. 

The fit is made using a maximum likelihood approach 
where the error bars are due to cosmic variance alone and 
are calculated from the model power spectra. The likelihood 
is 



C oc |Mr^''%xp{-(d^M"M)/2}, 



(13) 



assuming the covariance matrix M to be diagonal and 
therefore neglecti ng correlations between multipole bins 
(jBlake et al.ll200'i1 ). The vector d is of differences between 
model and data, di = C™""^"' - C^|*'' in each bin, and the 
diagonal elements of M are 



fsi.y{2£+l 



"C'e: 



(14) 



The fit is performed firstly for the photometric redshift data 
alone and secondly for the spectroscopic data alone. The 
best-fitting values for the bias are given in Tabled Both the 
spectroscopic and photometric data show an increase in bias 
with redshift, which gives us confidence that this effect is 
real. The difference in bias between all photometric redshifts 
and SDSS true-2 information is small in comparison with 
this trend; given the restricted overlap of these samples and 
the correlated errors shown in Fig. [5] we believe that it is 
insignificant within cosmic variance. 
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Galaxy angulor power spectrum for 0.0 < z < 0.1 



Galaxy angular power spectrum for 0.1 < z < 0.2 



Goloxy ongulor power spectrum for 0.2 < z < 0.3 






i 
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5 10 20 50 1 2 
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Figure 5. Fits to the measured galaxy angular power spectra in each redshift slice used to determine the bias. The bold symbols show the 
photometric data and the grey symbols the spectroscopic data (at ^ > 30 only) . Fits to the photometric data alone and the spectroscopic 
data alone are shown by the red solid and blue dotted lines respectively. 





0.0 < 2 < 0.1 


0.1 < 2 < 0.2 


0.2 < 2 < 0.3 


photo-2 only 


1.22 


1.65 


2.86 


SDSS only 


1.03 


1.90 


2.22 



Table 2. Best fit values for the bias calculated using the photo- 
metric redshift data only and the photometric data replaced by 
SDSS data using regions where the latter information is available. 



4 CROSS-CORRELATION ANALYSIS 

In this work, we have two hypotheses: a null hypothesis 
CgT = 0, and a hypothesis where the cross-correlation is 
as expected in a ACDM universe and is described by equa- 
tions (|8]) and (jlOp . We assume that the cosmological model 
is given; there are thus no adjustable parameters to consider, 
and the hypotheses can be compared by means of a simple 
likelihood ratio. 

The cross-correlation for a pair of galaxy and tempera- 
ture maps is calculated using 



CgT = 



1 



{21 + 1) 



(15) 



where the /sky factors attempt to compensate for the loss 
of sky coverage due to the masks. In detail, masking also 
introduces correlations between multipoles, which we take 
into account by using a full covariance matrix to calculate 



the X statistic; see Section [4.11 The cross-correlation data 
are binned to improve the signal-to-noise; we use 5 loga- 
rithmically spaced bins for 3 ^ ^ ^ 30 and ignore i = 2 
for consistency with R07. From our 3 redshift slices with 5 
multipole bins per slice, we have 15 data points in total to 
analyse for each WMAP band. 

The cross-correlations measured in each redshift slice 
for each of the WMAP bands are shown in Fig. (6] The re- 
sults are the same for each of the WMAP temperature maps 
considered, which is certainly consistent with the achromatic 
nature of the ISW effect. However, the plots also show the 



rms dispersion in CgT under the null hypothesis of no ISW 
effect, from which it is clear that no signal is significantly 
detected. We quantify this by evaluating the statistic for 
each hypothesis using 



(16) 



where C is the covariance matrix computed from simulations 
(see Section [4T]). The di are the values of (CgT " - C^^^) 
where is either zero or given by equations ([§} and H10|) 
for the fiducial hypothesis; CgT*^ are the measured cross- 
correlation values and i labels the binned measurements. 
The likelihood of each hypothesis is determined using 



£oc |Cr'/%xp|-(d'^C~M)/2|, (17) 

with d the data vector whose components are defined above 
and C the covariance matrix. In practice, the covariance ma- 
trix used in each case is calculated from simulations of the 
null hypothesis; this is a reasonable approximation given the 
subdominant nature of the ISW effect. Then 



21n^|i^=Ax^ 



(18) 



where £i and C2 are the likelihoods for the two hypotheses. 



4.1 Covariance matrix estimation 

In order to calculate the covariance matrix for the Cg-p val- 
ues, we use 500 simulations each of independent CMB and 
galaxy density fields. The power spectrum used for the CMB 
Gaussian realizations is the best-fitting theor etical angular 
powe r spectrum found by the WMAP team (jSpergel et al.l 
|2003). 

We generate simulations of the galaxy density fields by 
assu ming a lognormal form for the measured 2MASS density 
field (|Coles fc Joneslll99ll ). This allows a rapid generation of 
realistic density fields, from which mock galaxies are drawn 
via Poisson sampling. We now have a set of simulated CMB 
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0.0 < z < 0.1 0.1 < z < 0.2 0.2 < z < 0.3 




5 10 20 5 10 20 5 10 20 

1 I I 



Figure 6. The cross-correlation results for each of the photometric redshift slices using all the WMAP data: ILC (open triangle), Q 
(open circle), V (star) and W (filled triangle). We see from the Icr error bars around the null hypothesis (shown as lines) that the data 
are consistent with no ISW effect. 



0.0 < z < 0.1 0.1 < z < 0.2 0.2 < z < 0.3 




t I i 

Figure 7. The cross-correlation data for the V-band CMB data in the different redshift slices. Error bars are calculated from simulated 
CMB skies and 2MASS galaxy maps. The expected ISW signal in a ACDM universe is also shown calculated for bias values measured 
from the photometric data (see Table [2J . 



and 2MASS galaxy maps with the same shot noise properties 
as the data. 

The covariance matrix for the null hypothesis is esti- 
mated from the simulated data 

C,, = {(d, - {d.))(d, - (d,))) , (19) 

where the data vector d is defined as in equation (|16|l . Using 
only simulations of the CMB correlated with the observed 
galaxy maps leads to an underestimation of e rror bars at 
the 10 % level on all scales, in agreement with ICabre et al.l 
The error bars calculated from the diagonal elements 
of the covariance matrix are plotted in Fig.[7l and contrasted 
with the ISW signal expected in our default ACDM universe. 

4.2 Hypothesis testing 

The values for each hypothesis are given in Tables [S] H] 
and [S] considering bias values computed from photometric 
data only, SDSS data only and for b = 1.4 consistent with 
R07 respectively. Each table has data for each redshift slice 



and for an analysis of all three slices together. The results 
depend to some extent on the values adopted for the bias 
in each slice. This does not aff'ect the covariance matrix: 
values for the bias only affect the final results through 
an alteration of the expected CgT signal. 

In general, the results are low for the number of 
degrees of freedom; this was also noted by R07. We would 
expect values around = 5±vT0 for the individual redshift 
slices and x^ ~ 15±\/30 for all slices together, but only the 
second redshift slice gives x^ values that fall within these 
expected ranges. An analytic estimation of the error bars 
lAfshordi et al.ll2004) . 

'"'(^^t) = f d^U ^^'^-^ + QgCxT), (20) 

Jsky(^t. + ij 

suggests errors that are even larger than those used here 
except for the lowest I bin in the third redshift slice, so 
over-estimation of our error bars seems unlikely. 

For the total analysis using the biases calculated from 
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0.0 < 


z < 0.1 


0.1 < z 


< 0.2 


0.2 < 


z < 0.3 


0.0 < 2 < 0.3 


WMAP Band 


Xnull 


Xfld 


'^null 


Xfld 


Xiiull 


Xfld 


Xnull 




ILC 


0.70 


0.41 


3.8 


2.6 


0.52 


1.3 


5.0 


4.3 0.7 


Q 


0.68 


0.32 


3.7 


2.5 


0.68 


1.2 


5.1 


3.9 1.2 


V 


0.63 


0.33 


3.9 


2.8 


0.65 


1.4 


5.2 


4.4 0.8 


W 


0.63 


0.33 


3.6 


2.6 


0.67 


1.3 


4.9 


4.3 0.6 



Table 3. x^ values for the two hypotheses for each WMAP band analysed. The first three sections of the table show the values for 
the three redshift slices separately, the final section of the table shows x^ and Ax^ values for the entire data set. We see that the first two 
redshift slices prefer a ACDM ISW signal whilst the third prefers the null hypothesis. The whole dataset prefers a ACDM ISW signal. 
The significance of this result is discussed in 14.21 





0.0 < 


z < 0.1 


0.1 < z 


< 0.2 


0.2 < 


z < 0.3 


0.0 < 2 < 0.3 


WMAP Band 


Xnull 


Xfld 


^null 


Xfld 


Xiiull 


xid 


Xnull 


xid Ax^ 


ILC 


0.70 


0.44 


3.8 


2.6 


0.52 


0.86 


5.0 


3.9 1.1 


Q 


0.68 


0.36 


3.7 


2.4 


0.68 


0.82 


5.1 


3.5 1.6 


V 


0.63 


0.37 


3.9 


2.7 


0.65 


0.96 


5.2 


4.0 1.2 


W 


0.63 


0.37 


3.6 


2.6 


0.67 


0.96 


4.9 


3.9 1.0 



Table 4. As above except that bias values in the three redshift slices are calculated using spectroscopic measurements of the angular 
power spectrum (see Section 14.211 . 



(21) 



the photometric data alone, Ax^ = Xnuii ~ xid — 0.825, 
which translates to a likelihood ratio of 

£(lSW|data) ^ ^o.4l ^ ^ 
£(No ISW|data) ~ ^ - ■ ■ 

This is a long way from decisive evidence for the ISW 
effect. Using the other bias values we prefer the fiducial hy- 
pothesis at levels of 1.8 : 1 and 2 : 1 for the SDSS bias and 
b = 1.4 respectively. This latter scenario (6 = 1.4) is used 
by R07 and we find Ax^ values very similar to their results 
in this case: Ax^ — 1-5 against their Ax^ — 1-6. The differ- 
ences between the analyses are our use of photometric red- 
shifts, our slightly different magnitude cut of 12 < A" < 13.8 
rather than 12 < K < 14 and our use of simulations of the 
galaxy data as well as the CMB data to compute errors. 
As a result of simulating both CMB and galaxy data, we 
would expect our x values to be smaller, as indeed is the 
case. Thus the main result of this paper is that photometric 
redshift information does not increase the sensitivity of the 
2MASS data to the ISW effect. 



5 STATISTICAL POWER OF CURRENT AND 
FUTURE ISW EXPERIMENTS 

Given the lack of a significant ISW detection from 2MASS, 
one is led to ask whether this is as expected, or whether 
we have been unlucky. We can address this by using our 
simulation apparatus to generate the distribution of Ax^ 
values we could expect if the no-ISW null hypothesis were 
true 

P{Ax^ = x\Ho true). (22) 

Here, the statistic Ax^ itself is computed from the data in 
comparison with both Hq the null hypothesis, and Hi the al- 
ternative hypothesis but for data generated where Ho is true. 
This distribution therefore shows how frequently a given 



threshold in Ax^ would reject the null hypothesis when it is 
in fact true (the probability of making a Type I error). 

In a similar manner, we can obtain a distribution of Ax^ 
values under the assumption that the alternative hypothesis 
(ACDM) is true. For each pair of simulated CMB and galaxy 
maps, we add to the CMB map the expected ISW effect for 
the galaxy map in question. The expected ISW signal can 
be computed by using equation ([l| together with Poisson's 
equation to find: 



T 



= -2 



dg , 



a ^im[a) dr 
^3 



-Hia) ^(a) 



da 



n 



1 - d 



a 



'^im dr 



\-a)H{d)^-fAr, 



(23) 



l{l + l)c^ y "g(a) 

where 5 is the projected galaxy density field in the redshift 
shell under consideration and a = (1 -f- z)^^ , where z is the 
redshift at the midpoint of the shell. 

The set of simulated galaxy and CMB-I-ISW maps al- 
lows us to calculate the distribution of the Ay^ values we 
should expect if the alternative hypothesis were true: 

2 



P(Ax =x \ Hi true). 



(24) 



We can thus measure the probability of making a Type II 
error (accepting the null hypothesis when the alternative 
is true) for a given threshold. Fig. |8] shows the two Ax^ 
histograms computed from simulations of the CMB and 
2MASS galaxy data for the different redshift slices used in 
the analysis. 

A powerful statistical test would show a clear offset be- 
tween the Ax^ histograms in Fig. [8] In reality, we see very 
little offset between the Ax^ distributions for this exper- 
iment. The final panel of Fig. [5] reveals that if the ACDM 
hypothesis is true, a very large value of Ax^ would be needed 
to rule out the null hypothesis with any confidence. There 
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Table 5. As above but assuming a constant bias b = 1.4 in all three redshift slices as in R07. 
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Ax^ Ax^ 

Figure 8. The filled histograms show the values of Ax' calculated from simulation for the null hypothesis (no ISW effect); the outline 
histograms are for the alternative hypothesis (ACDM). Panels are 0.0 < z < 0.1 (top left), 0.1 < 2 < 0.2 (top right), 0.2 < 2 < 0.3 
(bottom left) and the full data set 0.0 < 2 < 0.3 (bottom right). A positive Ax^ indicates that the ACDM signal is preferred. Thus a 
powerful experiment would generate negative values of Ax' if the null hypothesis is simulated, and positive values if ACDM is simulated. 
In practice, we see little offset between the Ax^ distributions in all cases. As we move to higher redshift slices, the peaks of the Ax^ 
distributions begin to separate slightly, although any offset remains small. The broader histograms in the 0.2 < z < 0.3 slice are a result 
of shot noise in the galaxy density field. 



is an 8% chance of obtaining a Ax' value greater than 15 if 
the ACDM hypothesis is true, and a 5% chance of making 
a Type I error with this threshold. 

If we were to adopt Jeffreys' criterion (jjeffrevd Il948l ) 
that Ax' > 5 constitutes strong evidence for rejecting the 
null hypothesis, then such a threshold carries a ~ 23% 
chance of making a Type I error. 

To make a decisive detection of the ISW effect we need 
Ax' distributions with a clearer offset and less overlap. The 
poor 2MASS results prompt us to ask how good the data 



would have to be to make a significant ISW detection. As 
the CMB data on these scales are already signal-dominated, 
any improvement must come from the galaxy data. 
The power of a statistical test is defined as 

l~P{s.) (25) 

where P(s*) is the probability of making a Type II error with 
a threshold s*. A statistical test can be made arbitrarily 
powerful by altering the threshold, but a gain in power is 
offset by an increase in the probability of making a Type I 
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error. We therefore define the 'Optimal Power' of a test as 
the power of that test with threshold chosen such that the 
probabilities of making Type I and Type II errors are equal. 
A larger Optimal Power indicates a 'better' statistical test in 
the sense that it is less likely to return the wrong conclusion 
from the data. The Optimal Power of this experiment is 
0.55. 



5.1 Same survey, fainter magnitude limit 

Consider a survey with the same parameters as above but 
more galaxies at higher redshift (0.1 < z < 0.3), equivalent 
to having a fainter magnitude limit for the survey. Follow- 
ing the same procedure used in Section [4. II we simulate log- 
normal galaxy maps for this survey with bias and number 
density equal to the measured values for z < 0.1. 

The ACDM Ax^ distribution in the third redshift sUce 
(histograms not shown here) is narrower than that for 
the 2MASS experiment, although the 'null' distribution is 
largely unchanged, reflecting a reduction in shot noise which 
improves the ISW signal estimation. Looking at the analysis 
of the entire dataset together, the peak of the ACDM A^^ 
distribution is shifted to a higher value and the distribution 
is narrower, leading to an increase in the offset between the 
peaks of the distributions. However, the Optimal Power for 
this experiment is 0.60, showing little improvement. 



5.2 Deeper survey (more redshift slices) 

We also investigate hypothetical deeper surveys, with char- 
acteristic redshift Zm. The redshift distribution is computed 
using 



n{z) = ncom(z) exp{ — (z/zm)^'*"}. 



(26) 



where ncom is constant and calculated from 2MASS galaxies 
with z < 0.1. We also consider varying Zmax- Generally, we 
shall show results for a maximum redshift up to 0.7. 

We consider shells of width Az = 0.1 and take the true 
redshift distribution of galaxies in each slice as a Gaussian 
of width a — 0.03(1 + z) centred at the midpoint of the slice 
for z > 0.3; the true redshift distributions calculated for the 
2MASS data are used for z < 0.3. The galaxy bias in each 
slice is estimated using the simple assumption that bias is 
determined solely by the ratio of observed to expected num- 
bers of galaxies in that slice, where the expected number of 
galaxies is found using ncom. This means that bias is not 
intrinsically redshift dependent - the reason we see higher 
bias at higher redshift is because only the brightest, most 
massive galaxies are observed there. This seems a reason- 
able model for the r e dshift s we consider given the results of 
iMagliocchetti et al.l (|2000l ). who see little evolution in bias 
over this range. In keeping with these assumptions, we use 
the bias in the three 2MASS redshift slices, together with 
the fraction of galaxies observed here, to deduce a relation- 
ship between these quantities. Our best-fitting relationship 



b = 1.2/, 



gals 



(27) 



where /gais is the fraction of galaxies observed. 

Fig. |9] shows the Ax^ histograms for this experiment 
with Zm ~ 0.3. The Ax^ distributions for redshift slices with 



best possible experiment, z 



= 0.7 




0.5 



Figure 10. The build up of Optimal Power as a function of 
2max for the 2MASS experiment (red circles), more galaxies (blue 
stars), 2ni = 0.3 (green dashed line) and Zm = 0.75 (orange solid 
line). The top horizontal line shows the idealized case where the 
ISW signal (for Zmax = 0.7 is known exactly, rather than esti- 
mated from galaxy data. 



z > 0.3 all show an offset and the null hypothesis distribu- 
tions for these slices are slightly skewed towards negative 
values. The analysis of the full dataset for 0.0 < z < 0.7 
shows a clear offset between the peaks of the two distribu- 
tions and considerable improvement in detection prospects. 

Values of Ax'^ ^ 10 would be expected to occur (if the 
alternative hypothesis were true) around 30% of the time 
and would lead to a fairly conclusive rejection of the null 
hypothesis (5% chance of a Type I error). The Jeffreys cri- 
terion of Ax^ > 5 giving a 'strong' detection would have a 
9% chance of a Type I error and would be satisfied by the 
data 43% of the time under the ACDM hypothesis. The Op- 
timal power of this experiment is 0.73; Fig. [10] shows how 
this increases with depth. 

The results for an experiment with Zm = 0.7 are prac- 
tically indistinguishable from those for the 2m = 0.3 sur- 
vey and are therefore not shown in full. The only difference 
is a slight increase in offset in the Ax^ distributions for 
0.5 < z < 0.6 and 0.6 < z < 0.7. This is caused by a 
reduction in shot noise in these slices leading to a better 
estimation of the ISW signal. The Optimal Power of this 
survey is 0.74, its build up as successive redshift slices are 
added to the analysis is again shown in Fig. 1101 

The fact that the power depends mainly on Zmax and 
hardly on Zm may at first seem surprising. But as long as 
there are enough galaxies observed in a redshift slice that 
shot noise is small, and the bias in the slice is known, the 
number of galaxies and the particular value of b are unim- 
portant. 
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Figure 9. A simulated ISW detection experiment with characteristic depth Zm = 0.3 and Zmax = 0.7. The filled histograms show the 
values of Ax'^ calculated for the null hypothesis; the outline histograms are for the values of Ax'^ for the ACDM hypothesis. Comparison 
of the first three redshift slices with the previous 2MASS analysis shows improved detection prospects for these slices and overall. 
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5.3 Results for the 'ideal' case 

The best conceivable ISW experiment would measure di- 
rectly the temperature fluctuations due to the ISW effect 
(e.g. from perfect knowledge of the dark matter density field) 
across the entire sky, and cross-correlate such maps with the 
all-sky CMB. We consider this hypothetical limiting case 
in terms of our Ax^ analysis using realizations of the ISW 
angular power spectrum as our 'galaxy' maps. Even in this 
'best possible' case, we do not expect to be able to reject the 
null hypothesis with absolute certainty due to cosmic vari- 
ance on the scales of the ISW effect. The Ax'^ histograms 
we find for the ISW effect with Zmax = 0.7 are shown in 
Fig. [TT] and the probabilities of making Type I and II errors 
as a function of threshold in Fig. 1121 We notice that the 
separation of the distributions is fairly distinct, although a 
non-negligible overlap remains; the Optimal Power for this 
experiment is 0.95. 

Interestingly, if Ax^ ^ 5 constitutes 'strong' evidence 
for the ISW effect in a ACDM universe, then we would ex- 
pect to fail this criterion ~ 15% of the time i.e. only 85% 
of ACDM universes have an ISW effect that is detectable 
at this level. In the remaining ~ 15%, the intrinsic CMB 
signal serves to effectively 'hide' the ISW signal. Bear in 
mind that this experiment assumes all-sky measurement of 
the ISW temperature fluctuations to Zmax = 0.7, and that 
(necessary) masking of the sky will, as shown earlier, de- 
grade the power of the experiment with any realistic galaxy 
data, so that probably 1 in 4 observers would be incapable 
of detecting the ISW effect given data to Zmax = 0.7. 

This redshift limit of 0.7 includes about 50% of the ISW 
power compared to Zmax = oo. One can therefore do better 
with a larger redshift limit (e.g. Zmax = 1.3 would capture 
85% of the total power). But even in the limiting case of per- 
fect knowledge of the total ISW signal back to last scattering 
over the whole sky, we would fail to find 'strong' Ax^ 5S 5 
evidence for the ISW effect in 2% of ACDM universes. With 
the combined practical effects of realistic redshift coverage, 
galaxy sampling, and sky masking for the best foreseeable 
future surveys, this figure would rise to the order of 10%. 
Our results thus suggest that it is not unlikely that we could 
inhabit a universe where the ISW effect is present but un- 
detectable. 



6 CONCLUSIONS 

We have analysed photometric redshift data from 2MASS 
to Zmax = 0.3 and WMAP CMB data to look for a cross- 
correlation signal indicative of an ISW effect in a ACDM 
universe. This builds on the work of R07 who perform a 
similar analysis without redshift information. Our results 
are equally inconclusive: the data do not rule out a null 
hypothesis of no ISW effect; a ACDM ISW signal is pre- 
ferred, but with a likelihood ratio of only 1.5 : 1. We use 
error bars computed from simulations of both CMB skies 
and lognormal galaxy density fields and impose a slightly 
stricter magnitude cut than R07 on the 2MASS data to en- 
sure uniformity. The smaller sample size and larger error 
bars that we employ can explain the fact that our x^ values 
are slightly more unusually low than those of R07, but the 
generally low x^ nature of this dataset remains. 
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Figure 11. The filled histogram shows the Ax'^ values for the 
null hypothesis in the 'ideal' case where we could directly measure 
the ISW temperature fluctuations to Zmax = 0.7; the outline his- 
togram is for the ACDM hypothesis. We see a much diminished 
overlap between the distributions in this ideal case, although cos- 
mic variance means that there is not complete separation. 
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Figure 12. The probability of making a Type I error (solid line) 
and a Type II error (dotted line) as a function of threshold in 
Ax*^ for the ideal experiment with 2max = 0.7. The Optimal 
Power of this experiment (0.95) is 1 — p where p is the value of 
the probability where the above lines cross. We note that under 
the ACDM hypothesis, there is a ~ 15% chance of not detecting 
the ISW effect with a threshold for detection of Ax^ ?S 5. 



In the light of our inconclusive results, we have con- 
sidered the expected power of several ISW detection exper- 
iments. As it turns out, 2MASS is expected to be poor at 
discriminating between the two hypotheses considered. Sim- 
ulations of hypothetical deeper photometric redshift surveys 
reveal detection power to mainly depend upon the value of 
Zmax used for the analysis rather than the size of the galaxy 
sample. As long as sufficient galaxies are observed in each 
redshift slice to keep shot noise from dominating, the precise 
number of galaxies is largely irrelevant. 

The limiting case for such detection experiments is that 
of the 'ideal' ISW experiment where one can measure the 
ISW temperature fluctuations precisely out to a given Zmax. 
Even in this over-idealized case, we do not expect perfect 
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detection prospects due to cosmic variance on these scales: 
a certain fraction of intrinsic CMB skies have large-scale 
temperature fluctuations that 'hide' the ISW effect. Simu- 
lations indicate that, in this ideal case, a Jeffreys' criterion 
of Ax'^ = Xnuii ~ xid 5 for a 'strong' detection would fail 
to be satisfied ~ fO% of the time, for the best conceivable 
practical datasets. 

Given that the ISW effect seems doomed to be a low 
S/N phenomenon, is there a conflict with existing studies 
that claim to have detected the effect at high levels of signif- 
icance? Not necessarily; our results concern how frequently 
one will fail to find compelling evidence for the ISW effect, 
and do not say how strong the evidence will be if we are 
lucky enough to live in a universe where this detection is 
possible. Nevertheless, there is reason to doubt how much 
further the ISW effect can be pushed as a probe of preci- 
sion cosmology. In terms of sky coverage and redshift range, 
most of the ISW signal remains to be mapped. Even if we 
have been lucky with the regions of space studied to date, 
there is no guarantee that this will apply to the total sig- 
nal. Future studies of the ISW effect may fail to increase 
the significance of current detection s, and could even reduce 
them; the recent non-detection bv lSawangwit et al.l (|2010l ) 
using a high-redshift LRG sample may be an example of this 
phenomenon. 
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